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ABSTRACT 

We investigate numerically the effects of various boundary condi- 
tions on the flow field characteristics of the Physical Vapor Transport 
process. We use a prescribed temperature profile as boundary condition 
on the enclosure walls, and we consider parametric variations applicable 
to ground-based and space microgravity conditions. For ground-based 
applications, density gradients in the fluid phase generate buoyancy- 
driven convection which in turn disrupts the uniformity of the mass flux 
at the interface depending on the orientation. Heat conduction in the 
crystal can affect the fluid flow near the interface of the crystal. When 
considering isothermal source and sink at the interfaces, we observe a 
diffusive mode and three convective modes (i.e., thermal, solutal and 
thermo-solutal). The convective modes show opposing flow field trends 
between thermal and solutal convection; theoretically, these trends can 
be used to achieve a uniform mass flux near the crystal. However, under 
the physical conditions chosen, the mathematical condition necessary for 
uniform mass flux cannot be satisfied because of thermodynamic restric- 
tions. When a longitudinal thermal gradient is prescribed on the boun- 
dary of the crystal a non-uniform interface temperature results, which 
induces a symmetrical fluid flow near the interface for the vertical case. 
For space microgravity applications, we show that the flow field is dom- 
inated by the Stefan wind and a uniform mass flux results at the inter- 
face. 
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NOMENCLATURE 


a 

At 

b 

C P 

Dab 

f 


Constant term in the vapor pressure law (a=29.75) 
Aspect ratio, 


_L 

H 


Constant term in the vapor pressure law (b=11767.1K) 
Thermal heat capacity of the gas 
Binary diffusion coefficient 
Functional relationship 


F(x,y) Source and sink term. 


= jpC p u^dy 


g Gravitational acceleration 

go Gravitational acceleration at sea level 


Gr, Thermal Grashof number. 


gp 2 // 3 pA7 

P 2 


Gr, 

H 

Ja 


Solutal Grashof number, 
Height of the enclosure 


gp 2 H 3 aA(D A 


P* 


Dimensionless mass flux of A 




pu (L +1 ,y ) m „ 

i™, Dimensional mass flux of A ( mol/s-m 2> = pu(L+l,y) mtl ) 

Jr Total mass flux of A 


H 


J T = |p (L+l,y)u(L+l,y)dy 


k Thermal conductivity of the gas 

K Degree Kelvin 

L Vapor transport length 

/ Length of the crystal 
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Mi Molar mass of compound i 
P Pressure 

Pe Peclet number, In 

Pr Prandtl number, C P ^ 

< * 

R Perfect gas law constant (R = 8.32 kgm 2 s~ 2 mor i K~ l ) 

Sc Schmidt number, 

T Temperature 

u Velocity in the x-direction 

v Velocity in the y-direction 

x* ,y* Dimensionless coordinate axes 

Greek symbols 

a Coefficient of solutal expansion 

p Coefficient of thermal expansion 

\ Thermal conductivity of the solid 

p Dynamic viscosity 

co, Mass fraction of component i 

p Density of the gas 

0 Angle between the ampoule axis and the gravity vector 
Subscripts 

A Compound A 

B Compound B 

c Pertains to the crystal interface 

s Pertains to the source interface 

T Total 

x Component in x direction 

y Component in y direction 

Superscripts 

c Pertains to the crystal interface 

s Pertains to the source interface 

* Dimensionless 




>,(L+/)V 
Pad) . 
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L INTRODUCTION 

A sealed ampoule containing source material is typically used to grow large sin- 
gle crystals of opto-electronic materials by the process of Physical Vapor Transport 
(PVT). In the mathematical formulation of the PVT process a two dimensional enclo- 
sure has been considered For our physical situation, we also consider a 

rectangular enclosure, but we include conduction effects in the crystal. In this case, 
represented in Figure 1 , the nutrient is heated on one end, sublimes and transports 
toward the cooler end, where it condenses. This diffusive flow is known as advective 
flow or Stefan wind. However, since the process runs under an imposed thermal gra- 
dient, significant density gradients are possible; consequently, a naturally recirculating 
convective flow can occur in addition to this basic flow. These gradients can arise 
from both thermal gradients (causing thermal convection) and concentration gradients 
(inducing solutal convection). 

In recent years, publications on crystal growth from the vapor phase have mainly 
focused on the study of thermal convection 0-4,8,10,19^0) its on the non pianar 

interfacial mass flux <* ,1U,,20) . Most of the authors have included solutal convection in 
their formulation, with the use of the framework of the Boussinesq approximation; or, 
alternatively, prescribed their species boundary conditions. The study of the process 
with prescribed species boundary conditions allows the effect of thermal convection on 
the process to be studied independently < ^ 8 ' 18 “ 20) . Therefore, the predicted mass flux dis- 
tributions at the interfaces are influenced almost entirely by thermal convection. Rosen- 
berger et al. (4H> - U) investigated extensively combined thermo-solutal convection. But 
they were not concerned with the particular case of mercurous chloride for which solu- 
tal convection becomes dominant because of both high molecular mass and vapor pres- 
sure. Recently, Extremet et al. (18 > have studied simultaneously the balanced effects of 
thermal and solutal convection in a multizone rectangular enclosure. In their analysis, 
they used the Boussinesq approximation and they showed that special kinds of 
behavior, such as multiple counter-rotating cells, result from the opposing flow field 
trends between solutal and thermal convection. 

Heat transfer has a predominant role in convection occurring during crystal 
growth 09aa> . Previous studies have considered geometrically flat interfaces at a con- 
stant temperature. However, if there is a temperature gradient in the crystal, this could 
induce localized recirculating cells which would modify the mass flux along the inter- 
face. In addition, a curved crystal interface could itself play a role in inducing convec- 
tion cells. 

The main objective of this work is to study the effect of heat conduction in the 
crystal on the transport process during PVT for various imposed temperature profiles. 
We also investigate the opposing effects of thermal and solutal convection on the recir- 
culation patterns, in an attempt to exploit these contrary trends to correct nonuniformi- 
ties of mass flux at the growing interface in the case of combined thermo-solutal con- 
vection. In addition, we examine the mathematically necessary conditions which pro- 
vide for the existence of unidirectional flow in light of thermodynamic considerations. 
The microgravity environment is considered as a means of suppressing recirculating 
flow in an effort to restore the process to diffusion-controlled conditions. Finally, we 
investigate the role of longitudinal thermal gradients on the boundary of the crystal, 
especially as it relates to mass flux distribution at the interface and convection patterns. 


n. MATHEMATICAL MODEL 
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Following previous models, we consider a rectangular enclosure of height H, 
vapor transport length L, and crystal and source material length / , where / = L/2. This 
represents a two-dimensional cross-section of an infinitely wide PVT enclosure. The 
cavity’s axis is inclined at an angle 0 with respect to gravitational acceleration. The 
solid nutrient A exists at both sides of the ampoule in equal amounts. One side is at a 
temperature T s (source), the other at a temperature T c (sink, i.e., the growing inter- 
face). With T s > T c , the compound A sublimes from the source, is carried by a non- 
reacting buffer gas B, and condenses on the crystal interface. The source and sink 
interfaces are assumed to be geometrically flat. Their temperature is determined by a 
heat balance in the gas and the crystals. The temperatures of the source and the sink in 
contact with the enclosure walls and ends are either uniform in temperature or sub- 
jected to a linear thermal gradient A linear temperature profile from T s to T c is 
assumed along the upper and lower walls. The thermal conditions at the interfaces 
impose thermodynamically the weight fraction of component A in equilibrium with the 
solid through the partial pressure of the gaseous phase of compound A: 



We assume complete rejection of B on both interfaces, and bounding walls are 
impermeable to both components. The molecular weights of the two vapor components 
are different, specifically one being mercurous chloride (M A — All g/mol), the other 
being argon (M B = 40 g/mol). This will expand our study to include thermal as well as 
solutal convection. The viscosity, binary diffusion coefficient, and thermal conductivity 
are assumed to be concentration and temperature independent. 

The transport in this enclosure is governed by the system of nonlinear, coupled 
conservation equations for mass (continuity), species and energy, and balance equation 
for momentum (Navier-Stokes). These equations are, respectively, in rectangular two- 
dimensional coordinates, and for steady-state conditions: 

In the vapor, (/ < x <L+l ; 0 < y < H ) : 

Continuity equation: 

d(p«) + d(P v ) = 0 (2) 

dx dy 


Navier-Stokes equation 
Projection on the x axis: 


d . . 8 , . 3 P 

(pun) + — (pv«) = "fa 


dx 


dy 




a u 2 

du 3v 

dx 3 

a* + ay 

k / . 


a*v d 2 u 


(3) 


Projection on the y axis: 


d 


d , . a , . a p 

- (puw) + - ( pw) = -— 

a 2 " ^..3^ 


dy__2 

f ^ ’ 

du dv 

dy 3 

dx + dy 

k - 


+ Pgy 


(4) 
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Conservation of species equation: 


3(0, 

a 

3(0, 

pD "irj 


r* 


d , . 3 , . 3 

j£(pu <«,■)+ ^(pv»,)= — 


Conservation of energy in the fluid phase: 

|-(p»c,r) + | r (pvc,r)=|. 


and in the solid, (0 < x < l ; L+l < x < L+ 21 ; 0 < y < H) : 

3 2 r 3 2 r 
dx 2 by 2 


■* " 

k K 

.A 

\*L 

* dx i 

^ J 

3 y 



( 5 ) 


( 6 ) 


( 7 ) 


Boundary conditions: 

On the upper and lower walls, (/ <x <L+l, y=0 or y=H), the boundary conditions 
of impermeability, rigidity, no-slip and imposed thermal gradients were used: 


= 0 


3(o,(x,0) 3(0, (xfl) 

3 y 3 y 

v(x,0) = = 0 

u(x,0) = u(xjt) = 0 


3 T (x ,0) _ 37 (x,//) _ fs ~ Tc 
dx dx L 


( 8 ) 

( 9 ) 

00) 

(11) 


On the source, (x = l ; 0 < y < H), the diffusion of A from the source and impermeabil- 
ity to B were used, in addition to no-slip and an energy balance at the crystal inter- 
face: 


u(ly) = - 


E>ab 3(o a (l,y) 
l-(o A (/,>>) 3x 

v(/,y) = 0 




\dT(l,y) 


\ an/ ,y)' 

k dx 

V. J 

gas 

l J 


solid 


Analogously on the crystal, (x = L+l ; 0 <y <H) : 

Dab 3(O a (L+/0») 


u(L+l ,y) = — 


l-(0 A (L+l,y) 
v(L+l,y) = 0 


3x 


F(L+ly)~ 


' dT(L+ly) 


\d T(L+l,y) 

dx 

K J 

g<* 

K dx 

>> * 


solid 


( 12 ) 

03 ) 

(14) 

(15) 

(16) 
(17) 
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Boundary conditions for the 

(0 < x <, l ; L+l <, x 5 L+21 ; y = 0 or y = H) : 


source 


Case 1 


9r(*,0 ), = dnx jg l Ts-Tc_ Case2 


21 

T s -Tc 

l 


Case 3 


with corresponding temperatures at the ampoule ends: 

(x = 0 ; 0 S y ZH) , T = T, ; (x = 2L ; 0 <,y <>H) , T = T c 


crystal 


(18) 


(19) 


Most of these boundary conditions have been discussed in detail in previous pub- 
lications o/hu 1.12.14,18-2°) We w in therefore comment only the new ones we are taking 
into account. In particular, the matching conditions at the interface between the source 
material and the fluid phase yields equation 14 with a source term (F(l,y)) accounting 
for sublimation; likewise a sink term (F(L+l,y)) also appears in equation 17 for the 
condensing interface at the crystal end. Note, in order to satisfy continuity the source 
and sink terms are compatible. Since the growth rate of the interface is much smaller 
than the vapor velocities, typically by three orders of magnitude, the process is quasi- 
steady. Thus, the term accounting for latent heat of sublimation is negligible incom- 
parison to the conduction and advection terms. Also, similar to previous investigators 
we are considering that the interface remains flat during the growth of the crystal. 


ffl. NUMERICAL FORMULATION 

The system of coupled nonlinear partial differential equations has been solved 
using the SIMPLE algorithm (Semi-Implicit Method for Pressure-Linked Equations 
(22) ). The conservation equations were discretized using a hybrid scheme, which pro- 
vides a second-order spatial accuracy. This two-dimensional version of the algorithm 
uses three staggered grids, one for the scalar variables and one for each of the velocity 
components. In contrast to a vorticity/stream function approach, adopting a primitive 
variable approach allows a straightforward formulation of the velocity boundary condi- 
tions. However, the solution of the continuity equation for this problem is not automat- 
ically satisfied. 

For most of the simulations, a 20x40 primary grid was used. An additional case 
was run with a fine 40x60 mesh. The code was run on VAX 8800 biprocessor com- 
puter, and the iterative process was stopped when a Cauchy-Schwartz type criterion 
was satisfied for each variable. All the results presented in this article have been 
obtained with a convergence criterion of KT 4 for each of the variables. The code was 
compared with a more common vorticity/stream function algorithm under the same 
conditions for various Grashof numbers. Results provided by the two codes differed by 
less than one percent (21) . 

IV. RESULTS AND DISCUSSION 

Due to the finite value of crystal thermal conductivity, heat conduction in the cry- 
stal may affect the interface temperature and local convective patterns. We consider 
two cases: specifically, the flow field behavior when heat conduction is absent, and the 
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effect of heat conduction when longitudinal thermal gradients are applied at the crystal 
boundary. In the absence of heat conduction, we study the effects of the aspect ratio 
and the orientation of the ampoule with respect to the gravitational acceleration on 
thermal convection. Thermal, solutal, and thermo- solutal convection are studied in the 
case of a horizontal enclosure. We conclude this first part with a feasibility study of a 
unidirectional flow. Then, we show in the second case that when heat conduction in 
the crystal is present the local convection patterns are modified. However, under 
microgravity conditions, even in the presence of heat conduction, the flow field is 
restored to a unidirectional pattern. 

The governing parameters of this process, i.e., the thermal and solutal Grashof 
numbers, can be obtained by linearizing the density field using the Boussinesq approxi- 
mation: 


P = Po (1 - a(( 0 A -tfl£) - p (T-Td) 
where the subscript refers to these conditions 


T 0 = 


T c +T S . A _ cof + o oj 


and Wq = 


( 20 ) 


( 21 ) 


The coefficients of solutal and thermal expansion are derived from the ideal gas 
law and may be written as follows: 




2 l 

Tc + T s T 0 


( 22 ) 



P 


8p _ M s -M a 

J P T 


(23) 


Typical operating conditions for crystal growth of mercurous chloride by Physical 
Vapor Transport (PVT) are shown in table 1. The vapor pressure, density and mass 
fraction are calculated from equations 1, 25 and 30 respectively. The total operating 
pressure is taken as P to , = 20000 Pa; the excess pressure takes into account the effect 
of the residual buffer gas. 



1 









The dimensionless parameters for the process under the conditions of this study are 
shown in table 2. 
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Table 2: Typical dimensionless parameters | 

Thermal Grashof number = 

2700 

Solutal Grashof number = 

10000 

Prandtl number = 

0.7 

Schmidt number = 

1 

Aspect ratio = 

10 

1 < Mass Peclet number < 30 | 


cTT 

IV. 1. Constant temperature imposed on the crystal (Case 1, — = 0 ) 

IV. 1.1. Effect of aspect ratio on interface temperature 

We examine the effect of the aspect ratio on the interface temperature. The 
imposed thermal profile along the boundary of the enclosure is shown in Figure 2 
(case 1). Parametric studies for various aspect ratios are shown in Figure 3. These 
cases are computed for the horizontal configuration where convection is at a max- 
imum. The results show that in both the source material and the growing crystal, the 
interfacial temperature deviates more severely from a constant value as the aspect ratio 
decreases. The maximum temperature deviation (1.25 K) occurs for a square ampoule 
(Ar=l). However, for narrow ampoules with aspect ratio > 10, the interfacial tempera- 
ture approaches a constant value. 


IV. 1.2. Thermal convection in a rectangular inclined enclosure 

We now investigate the flow field characteristics occurring during PVT. We sim- 
plify the calculation by assuming isothermal source and sink interfaces. This 
corresponds to the absence of heat conduction in the crystal. We first study exten- 
sively thermal convection. 


IV. 1.2.1. General pattern of thermal convection 

Solutal convection is eliminated mathematically by setting the molar masses of 
both components to be equal: this implies ot = 0. Once we have set M A -Mb> we can 
use a simplified formula for the density given by the ideal gas law: 


M P 
p_ R T 


(24) 


Using this formula to compute density of vapor on both interfaces gives: 
p 5 = 1.93 kg/m 2 and pc =2.11 kg/m 3 . 

The intensity of thermal convection can be quantified using the thermal Grashof 
number defined as follows: 

G is^m. 

P 

The pattern of convection for Gr, = 2672 is shown in Figure 4.a, which exhibits 
the well-documented result of purely thermally driven fluid flow. The gas, lighter on 
the hot side, rises up along the nutrient interface, and heavier on the cold side, moves 
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down along the interface. Notice that the mass flux distribution. Figure 4.b, on the cry- 
stal interface exhibits significant nonuniformity. As is typical for thermal convection, 
more nutrient ascends to the top of the ampoule than is transported to the bottom of 
the ampoule. Figure 4.c shows that the convection stream modifies the mass distribu- 
tion throughout the entire ampoule. Note that when M A =M B , © ranges from .17 to .72. 
We conclude that thermally driven convection is the leading cause of nonuniform mass 
flux distribution for the case considered. 

IV. 1.2.2. Effect of the aspect ratio on thermal convection 

From a crystal growth point of view, we would like to reduce the convection 
intensity in order to obtain a uniform mass flux next to the interface. One way to 
reduce the effect of convection is to modify the geometry of the ampoule. Aspect 
ratios varying from 0.5 to 20 have been studied. Figure 5 presents the velocity fields 
resulting from such simulations. It demonstrates that the higher the aspect ratio (length 
to width ratio), the lower the convection intensity. With decreasing aspect ratio, 
viscous effects become important and a transition to a uniform flow field results. The 
transition from a convective pattern (with a recirculating cell) to a diffusive pattern 
(uniform flow field without a recirculating cell), is observed when the aspect ratio 
changes from 10 to 20. The variation of mass flux with aspect ratio is shown in Figure 
6. The maximum occurs for Ar=2, beyond which there is a sharp drop, and finally lev- 
els off. 


IV. 1.2.3. Effect of the angle between the ampoule axis and the gravity vector 

The angle between the ampoule and the gravitational acceleration vector is 
another control parameter of the process. We have systematically studied its effect for 
angles varying between 0 0 (thermally stable configuration) and 180 0 (thermally 
unstable configuration). We operated below the critical Grashof number. As a result, 
these trends should be viewed as only qualitative trends, in particular, the vertical 
configurations. 

The resulting velocity fields are represented in Figure 7. Both vertical 
configurations correspond to a diffusive mode of mass transport, since we are operat- 
ing below the critical Grashof number. As soon as the ampoule is tilted, by as little as 
one degree, a recirculating cell appears. The intensity of the convection increases with 
slight angular variation, and peaks for an angle of 90 °. It should be emphasized that 
we are obtaining qualitative data from this simulation because of possible three- 
dimensional effects. 

The distribution of the total mass flux for various angles of tilt is represented in 
Figure 8. The total mass flux variee between a minimum quantity, associated with the 
vertical configuration, and a maximum value for the horizontal configuration, which 
corresponds to the peak in convection levels. The minimum of the curve is sharper 
than the maximum, denoting the immediate transition to convection upon inclining the 
ampoule. Note that the total mass flux differs approximately in a proportion of 1 to 2 
between the diffusive and convective modes. The mass flux distribution along the 
interface for various angles of tilt is shown in Figure 9. Only the vertical configuration 
(0 °) of the ampoule results in a uniform mass flux. This is consequently the only case 
in which planar growth could possibly occur. 
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IV. 1.3. Solutal convection in a horizontal enclosure 

We have already noticed from the magnitude of the thermal and solutal Grashof 
numbers (table 2) that solutal convection dominates thermal convection. However, to 
isolate the effects of solutal convection, the temperature difference between the source 
and the sink is set to zero, keeping the mass fraction difference equal to what was pre- 
viously used. This is of course unrealistic from a physical standpoint, but our goal is to 
study separately each of the basic phenomena involved in crystal growth. 

The formula of the density derived from the perfect gas law cannot be simplified 
further, and the molar mass has to be written explicitly: 

i 

P 

( 25 ) 

The computed densities on the interfaces are: 

p, = 1.68 kg/m 3 and p c = 0.41 kg/m 3 

The intensity of solutal convection can be quantified using the solutal Grashof number 
defined as follows: 


P = 


m A 


(l-co A ) 


m B 


R T 


gp 2 H 3 aA( 0 A 


Unlike thermally driven convection, the density is greater at the hot interface than 
at the cold interface. We therefore expect a fluid flow and a mass flux distribution 
opposite to that observed for a thermally driven fluid flow. Figure lO.a presents the 
flow field resulting from such a theoretical hypothesis. The fluid flow moves down 
along the hot interface, and up along the cold interface. The maximum velocity is 
larger than that observed by purely thermally driven flow, indicating a stronger recircu- 
lation. This can be attributed to the larger magnitude of the Grashof number (Gr, of 
order 10 s ). The interfacial mass flux distribution shows that more nutrient flows to the 
bottom of the crystal interface than to the top of the interface (Figure lO.b), where the 
convective flow opposes the advection. 


IV. 1.4. Thermo-solutal convection in a horizontal enclosure 

Physically, one cannot of course set thermal or solutal convection to zero; both 
are present simultaneously. With this cautionary note in mind, we proceed to choose 
the temperature difference between the interfaces to be 50 K, which is representative 
of an average growth condition. Computing densities on both interfaces gives: 

p s = 1.61 kg/m 3 and pc = 0.42 kg/m 3 

The solutal Grashof number is approximately four times greater than the thermal 
Grashof number (9856 vs. 2672, respectively). Thus the flow field is of solutal convec- 
tion type. The mass flux distribution on the growing interface is therefore expected to 
be characteristic of solutally driven growth. 

We have shown that, in the case we are modeling, two types of convection are 
competing. Our goal being to restore a uniform mass flux, we can exploit these oppos- 
ing trends to eliminate their effects in the neighborhood of the growing interface. This 
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could theoretically be achieved if we set the Gr, equal to the Gr,, leading to the simple 
condition: 


P AT = a Aa A (26) 

This can be realized by increasing the temperature difference and decreasing the mass 
fraction difference between the two interfaces. However, this condition is not found to 
be sufficient to restore a uniform mass flux. A modification to equation (26) is neces- 
sary to account for the difference in mass-flux and thermal gradients; the thermal gra- 
dient of the furnace is imposed to be linear, but the mass-flux gradient has no such 
restriction. Most of the mass-flux gradient is localized near the interfaces, especially 
near the growing interface, so that the following arbitrary criterion becomes necessary: 

P AT = 2 a Aco* (27) 

The resulting velocity field is presented in Figure 11. a. A strong thermal recircu- 
lating cell is present throughout most of the ampoule, with a weaker recirculating 
region near the interface. Higher resolution studies (Figure ll.c) show that this 
weaker cell is unidirectional. The unidirectional flow does not occur throughout the 
ampoule, but is apparent close to the interface. This is enough to enable the Stefan 
wind to be the locally dominant mode of mass transfer. This results in a nearly planar 
mass flux distribution along the growing interface as shown on Figure ll.b. 


IV. 1.5. Feasibility of a unidirectional flow. 

To achieve a flow by satisfying the pseudoequality of the thermal and solutal 
Grashof numbers, we had to decrease the difference in mass fraction and increase the 
difference in temperature between the two interfaces independently of the thermo- 
dynamical laws regulating the process. The mass fraction on the interfaces is uniquely 
a function of the temperatures of the interfaces through the vapor pressure law of the 
solid/vapor equilibrium between the solid and gas phases of component A. The 
pseudo-equation of equality of the thermal and solutal Grashof numbers should then be 
written: 


p AT = 2 a A(o a 
m-T c ) - 2 a (coi-atf) 


(28) 


The thermodynamical vapor pressure law (equation 1) is utilized to produce the 
mass fraction through the following conversion equation: 


(0 A = 


PaMa 

PaM a + ( Pt-PaWb 


which can be rearranged to 


g>* = 



(29) 


(30) 


substituting equation (1), yields: 
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< »A = 


1 


(31) 


1 + 


P T e 


-1 


M b 

Ma 


If we substitute this formula in equation (28), we obtain: 


P (Ts -T c ) - 2a/ ,(r s Jc) = o (32) 

In other words: 

f 2 (T s J c fMAMB>Q (33) 

We will apply this functional relationship to our particular case with the following 
molecular masses of HgCl 2 and Ar, and total pressure in the ampoule: 

Table 3 1 

P = 20000 Pa 
M a = 0.472 kg/mol 
Mb = 0.04 kg/mol 

We also impose the following conditions on T s and T c : T S =T 0 +AT , T C =T 0 -AT with 
7 * 0 = 562 /:, corresponding to the working temperature of the process. The following 
simplified equation has to be solved: 

/ 3 (AI) = 0 (34) 

The solution shows that only one trivial root exists for AT=Q. For this root, there is no 
driving force for the process. We conclude that, under the conditions for which we stu- 
died the process, a unidirectional fluid flow using competition between thermal and 
solutal convection cannot be achieved in a terrestrial environment because of a thermo- 
dynamical restriction. Other working conditions for the process have to be found. 


iV.2. Imposed longitudinal thermal gradient in the crystal (Cases 2 and 


dT 

dx 


> 0 ) 


3, 


IV.2.1. Heat conduction in the crystal, 


dr 

dx 


>0 


The effect of heat conduction in the presence of thermal gradient on the crystal 
boundary (cases 2 and 3 in Figure 2) is shown in Figure 12 for the vertical 
configuration. These trends show that the steeper gradients increase the distortion in 
interface temperature. The corresponding flow field configuration for the steepest gra- 
dient (case 3) is shown in Figure 13.a. For the set of parameters chosen, the tempera- 
ture gradients at the interface give rise to symmetrical cells near the interface. The 
core flow is the most intensive, which results in a bell shaped distribution of mass flux 
(Figure 13.b). In contrast to our results for the horizontal configuration which showed 
slight temperature deviation from a constant value at the crystal interface (case 1 ), the 
vertical configuration shows that the interface is isothermal. These slight deviations can 
be accounted for by the large difference in Grashof number between the two 
configurations. In the calculations for the vertical configuration, we selected very low 
temperature differences (A 7 '= 1 . 5 A') to examine a limiting case for the PVT process. 
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IV.2.2. Microgravity environment 

One of the possible alternatives is to grow the crystal in a microgravity environ- 
ment. We have simulated different levels of gravity. The resulting fluid flow is shown 
in Figure 14 for KT 3 g 0 (case 3). For gravity levels of less than KT 3 ^ 0 (i.e., Gr < 10), no 
recirculating cell occurs. The corresponding mass flux distribution in front of the grow- 
ing interface is therefore planar. 

V. CONCLUSIONS 

Considerable progress has been achieved in the past years in modeling PVT 
processes. Our study covers the flow field characteristics during the growth of Hg 2 Cl 2 
crystals in a rectangular ampoule under terrestrial and pg conditions for various 
imposed thermal gradients. Taking into account the real thermodynamical dependency 
of the mass fraction as a function of temperature reduces considerably the control we 
have on the process. Solutal convection becomes the main phenomenon governing the 
fluid flow inside the ampoule. In order to achieve our goal of a uniform mass flux near 
the interface (planar growth), we decided to study systematically each of the 
phenomenon involved in the process. 

We showed that: 

1. In the absence of conduction in the crystal, the interface temperature becomes 
uniform for Ar > 10. Otherwise, in the presence of conduction when linear gra- 
dients are imposed on the crystal boundary, the interface temperature becomes 
non-uniform, which gives rise to multicellular convection in the flow field. The 
temperature distortion is proportional to the thermal gradient imposed on the cry- 
stal. 

The flow field characteristics during PVT were studied by isolating the trends of 
thermal, solutal convection, and the combination of thermo-solutal convection. We stu- 
died thermal convection by artificially eliminating solutal convection, i.e., by setting 
the masses of both gas components equal. We have also shown that: 

2. The magnitude of convection decreases with the aspect ratio (Ar). Up to an Ar of 
10, a recirculating cell is observed. For an Ar of 20, the fluid flow switches to a 
diffusive mode, associated with a restored planar growth. 

3. The magnitude of convection varies approximately with the sine of the angle 
between the ampoule and the gravitational acceleration, as does the mass flux dis- 
tribution on the growing interface. The maximum mass flux occurs for the hor- 
izontal configuration. Only vertical configurations (below the critical Grashof 
number) provide potential for a diffusive mode of transport of mass, resulting in a 
uniform mass flux near the interface (planar growth). 

We then studied solutal convection by restoring the mass difference and setting 
the temperature difference between the interfaces to zero. We have observed that for a 
horizontal configuration: 

4. The general behavior of the fluid flow is reversed. The sense of recirculation of 
the flow is opposite to the one observed for thermal convection. Likewise the 
mass flux distribution on the growing interface follows a similar trend, where 
more nutrient comes on the bottom of the ampoule. 

We then studied the more realistic case of a combined thermo-solutal convection 
in the horizontal enclosure, finding that: 
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5. Solutal convection was the dominant mode of the fluid flow in the ampoule. 

The opposing trends of thermal and solutal convection led to the investigation of 
the possibility of obtaining a unidirectional flow. This was achieved by increasing the 
temperature difference between the two interfaces and simultaneously decreasing the 
mass fraction difference, independently of the requirements of the vapor pressure law. 
This analysis indicated that: 

6. A unidirectional flow near the crystal interface is possible when the thermal 
Grashof number is equal to twice the solutal Grashof number under the operating 
conditions chosen. 

We looked for the conditions of such a unidirectional flow within the framework 
of the thermodynamic requirements: 

7. For the pressure, molar masses and temperature we are studying, such a unidirec- 
tional flow is impossible. Solutal convection always overwhelms thermal convec- 
tion because of the disparity in Grashof numbers. 

8. Microgravity appears to be an efficient way of resolving the problems. For grav- 
ity levels of less than 10" 3 « 0 , the Stefan wind drives the flow. No recirculating cell 
can be observed, and a uniform mass flux is obtained near the interface. 

Buoyancy-driven convection in PVT is almost unavoidable at 1-g because of the 
density gradients which are omnipresent. What are the possible actions at 1-g? Two 
actions come to our mind: 

1. Choice of the buffer gas: most of our problems arise from the huge difference in 
molar mass between argon and mercurous chloride (1 to 12). We could use a 
heavier gas. 

2. We can change the geometry of the ampoule, using the conclusion of the study of 
the effect of the aspect ratio (narrow ampoule). 

Alternatively, the following may also be possible: 

3. We can modify the partial pressures of the components. Increasing the partial 
pressure of component B would reduce solutal convection, but also reduce the 
rate of growth, the diffusion coefficient decreasing with pressure. 
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Figure 1— Physical vapor transport enclosure. 



Figure 2.— Temperature prolile along the ampoule. 
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Figure 5.— Eftect ot the aspect ratio on the fluid flow (For Ar = 1 , 2, 5, 1 0, and 20 the Figure 6.-Mass flux versus aspect ratio (J max - .25x1 0 1 mol/s-m 2 ) 
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3.34 10 2 respectively). 
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Figure 8. — Effect o( the angle of tilt on the total mass flux inside the 
enclosure (Ar * 10). 
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Figure 9 — Effect of the angle of tilt on the nonuniformily of the mass flux 
on the growing interface (Ar « 10, J max - .38x1 O* 1 mol/s m 2 ). 



(a) Fluid flow dnven buy solutal convection (Gr t = 0, Gr s * 9856. Gr^G^ * «, Ar = 10). 
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(b) Corresponding mass flux distribution along the growing interface (J ma „ » 1 .5x1 O'* 
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(b) Corresponding mass flux distnbution along the growing interface 
(J max = 1.35x1 O’ 1 mol/s-m 2 ). 

Figure 11. 



FIGURE lie. 
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Figure 12.— Resulting temperature profiles of the crystal interlace for various 
temperature gradients, cases (1), (2), and (3) correspond to 2ll * 0, -2.5, and 
-5 0 Kcm‘ 1 respectively. 


22 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden tor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLY ( Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

December 1991 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Thermal and Solutal Convection With Conduction Effects Inside a 
Rectangular Enclosure 


6. AUTHOR(S) 


Christophe Mennetrier and Walter M B. Duval 


WU- 674 - 21 - 05 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E -5521 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105371 


11. SUPPLEMENTARY NOTES 

Christophe Mennetrier, Visiting Research Associate, Ecole Nationale Superieure des Mines, France at NASA Lewis 
Research Center. Walter M B. Duval, NASA Lewis Research Center. Responsible person, Walter M B Duval 
(216)433-5023. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Category 34 


13. ABSTRACT (Maximum 200 words ) 

We investigate numerically the effects of various boundary conditions on the flow field characteristics of the Physical 
Vapor Transport process. We use a prescribed temperature profile as boundary condition on the enclosure walls, and 
we consider parametric variations applicable to ground-based and space microgravity conditions. For ground-based 
applications, density gradients in the fluid phase generate buoyancy-driven convection which in turn disrupts the 
uniformity of the mass flux at the interface depending on the orientation. Heat conduction in the crystal can affect the 
fluid flow near the interface of the crystal. When considering isothermal source and sink at the interfaces, we observe a 
diffusive mode and three convective modes (i.e., thermal, solutal and thermo-solutal). The convective modes show 
opposing flow field trends between thermal and solutal convection; theoretically, these trends can be used to achieve a 
uniform mass flux near the crystal. However, under the physical conditions chosen, the mathematical condition 
necessary for uniform mass flux cannot be satisfied because of thermodynamic restrictions. When a longitudinal 
thermal gradient is prescribed on the boundary of the crystal a non-uniform interface temperature results, which induces 
a symmetrical fluid flow near the interface for the vertical case. For space microgravity applications, we show that the 
flow field is dominated by the Stefan wind and a uniform mass flux results at the interface. 


14. SUBJECT TERMS 


f 

Physical vapor transport; Crystal growth; Thermo-solutal convection 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 


NSN 7540-01-280-5500 


15. NUMBER OF PAGES 

24 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 













